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Multidrug resistance (MDR) contributes to the failure of chemotherapy and high mortality in non- small cell 
lung cancer (NSCLC). We aim to identify MDR genes that predict tumor response to chemotherapy. 199 
NSCLC fresh tissue samples were tested for chemosensitivity by MTT assay. cDNA microarray was done 
with 5 samples with highest resistance and 6 samples with highest sensitivity. Expression of ABCC3 mRNA 
and protein was detected by real-time PCR and immunohistochemisty, respectively. The association 
between gene expression and overall survival (OS) was examined using Cox proportional hazard regression. 
44 genes were upregulated and 168 downregulated in the chemotherapy-resistant group. ABCC3 was one of 
the most up-regulated genes in the resistant group. ABCC3-positive expression correlated with lymph node 
involvement, advanced TNM stage, more malignant histological type, multiple-resistance to anti-cancer 
drugs, and reduced OS. ABCC3 expression may serve as a marker for MDR and predictor for poor clinical 
outcome of NSCLC. 

Lung cancer is the leading cause of cancer morbidity and mortality worldwide. Non-small cell lung cancer 
(NSCLC) accounts for 80-85% of lung cancer with a 5-year survival rate of 15%\ Although surgery is 
regarded as the best possible treatment for NSCLC, only 20-25% of tumors are suitable for potentially 
curative resection 2 and despite complete surgical resection, the recurrence rate remains high (30-70%). 
Adjuvant chemotherapy may reduce recurrences by eradicating the subclinical and micro -metastatic disease 2 . 
However, the heterogeneity of NSCLC tumors and resistance to multiple drugs results in a poor response to 
chemotherapy and reduced survival 2 . Multidrug resistance (MDR) is a multifactorial process and a major obstacle 
to treatment for NSCLC patients 2 . Identification of biomarkers for predicting an individual's resistance to 
chemotherapy and for rationally selecting the most efficacious therapeutic agents is dependent upon under- 
standing the molecular alterations that contribute to MDR in lung cancer. 

In vitro anticancer drug sensitivity tests can identify biomarkers to predict chemosensitivity, improve che- 
motherapeutic response and design individualized chemotherapy 3 " 5 . Such tests as MTT (3-(4,5-Dimethylthiazol- 
2-yl)-2,5-diphenyltetrazolium bromide), CD-DST (collagen gel droplet- embedded culture drug sensitivity test) 
and HDRA (histoculture drug response assay) have been developed for a variety of malignancies. MTT chemo- 
sensitivity test might help to predict tumor response in ovarian cancer and identify the optimum adjuvant 
chemotherapy for esophageal squamous cell carcinoma. HDRA might be feasible for predicting the efficacy 
and selecting the proper anticancer drug for lung cancer 6 . CD-DST test for "new generation" anticancer drugs 
by using surgically resected specimens of NSCLC patients may predict tumor response to postoperative chemo- 
therapy 7 . We applied MTT assay using resected NSCLC tumor tissues to stratify patients into chemotherapy- 
sensitive and resistant groups. We then analyzed the gene expression profile of the resistant and sensitive group by 
cDNA Microarray to identify genes associated with MDR. 212 genes with differential expression were identified 
including 44 upregulated and 168 down regulated in the chemotherapy-resistant group compared to the sensitive 
group. ABCC3 was one of the most up-regulated genes in the resistant group. 

ABCC3 is a member of the ATP-binding cassette (ABC) transporters family. ABC transporters are highly 
expressed in tumor cells where they actively efflux a broad spectrum of anticancer drugs and thus contribute 
to MDR 811 . ABCC3 is a member of the multidrug resistance-associated protein (MRP) subfamily which is 
involved in MDR. ABCC3 expression has been shown to be higher in NSCLC than that in SCLC 12 and 
ABCC3 expression correlated with decreased sensitivity of lung cancer cells to anticancer drugs (vincristine, 
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etoposide, and cisplatin) 13 , and especially to methotrexate and dox- 
orubicin in NSCLC 14 ' 15 . ABCC3 is also overexpressed in Anip973/ 
NVB cells consistent with involvement of ABCC3 in NVB-resistance 
in lung cancer 16 . However, the mechanism underlying the role of 
ABCC3 in MDR to taxanes (paclitaxel and docetaxel), gemcitabine, 
and cisplatin in NSCLC remains largely unknown. To assess whether 
ABCC3 expression has predictive value in chemosensitivity and pro- 
gnosis in NSCLC, we examined the expression of ABCC3 in patients 
with resected NSCLC tumors and evaluated its association with 
tumor characteristics and sensitivity to chemotherapy. Our goal 
was to examine the specific genetic traits of tumors as a means of 
predicting their chemosensitivity to "new generation" anticancer 
drugs. The results, if confirmed, would possibly identify biomarkers 
useful to predict clinical outcome of NSCLC. 

Results 

MTT assay. MTT assay was successfully performed using primary 
cancer cells from tissues of 199 patients; 69 were highly sensitive to all 
the 5 drugs and were categorized as the highly sensitive group; 44 
were resistant to all the 5 drugs and were categorized as resistant 
group; and 86 were partially sensitive to anticancer drugs and 
categorized as the moderate-sensitive group. Samples for 11 
patients (6 patients with highest sensitivity and 5 patients with the 
highest resistance to all the 5 drugs) were selected for gene expression 
microarray. Resistance index of the 11 patients is shown in 
Supplemental Table 1. 

Gene expression microarray. To identify potential predictors for 
chemosensitivity, cDNA microarray was used to analyze the gene 
expression profiling of the resistant and sensitive group. Fold change 
measurements were analyzed by SAM. The transcripts with fold 
change >2 (upregulated) or <l/2 (downregulated) and |Score(d)| 
> 2 were considered to have significantly differential expression 
between the two groups. The differentially expressed elements 
were identified by the two-sided statistical tolerance interval (95%) 
with a cut-off of false discovery rate at 10%. Two hundred-twelve 




Figure 1 | Heat map generated by hierarchical clustering of differentially 
expressed genes between the 5 resistant and 6 sensitive samples (Rows: 
tumor samples; Columns: genes). Differentially expressed genes (Red: 
upregulated; Green: downregulated) were selected at false discovery rate of 
10%. 



(212) genes were differentially expressed between the resistant and 
sensitive group. Compared to sensitive group, 44 were upregulated 
and 168 downregulated in chemotherapy-resistant group. Those 
genes spanned distinct functional families involved in biological 
processes including cell apoptosis, cell cycle, cell adhesion, signal 
transduction and ABC transporter (Supplemental Table 2 and 3). 
Hierarchical clustering was performed on the 1 1 arrays based on the 
differentially expressed genes (Fig. 1). 

ABCC3 expression by real time-PCR in NSCLC. The upregulated 
genes with potential predictive value for chemosensitivity in NSCLC 
were evaluated by real time-PCR (date was not shown). ABCC3 was 
one of the most differentially expressed in the resistant group. We 
performed quantitative RT-PCR for ABCC3 in the 199 NSCLS 
tumor tissues. ABCC3 mRNA level was higher in the resistant 
group compared to the sensitive group (mean ± SD: 0.4587 ± 
0.1324 vs. 0.1324 ± 0.0866, P = 0.001) (Fig. 2). 

ABCC3 expression by IHC in NSCLC and its correlation with 
clinicopathological characteristics ABCC3 immunoreactivity was 
present in the cell membrane and cytoplasm of NSCLC tumor 
(Fig. 3). ABCC3 -positive expression was observed in 63 cases includ- 
ing 23 with squamous cell carcinoma (SCC), 38 with adenocarci- 
noma (ADC), and 2 with other types (Table 1). ABCC3 expression 
correlated with pN status (P < 0.001), pathological TNM stage (P < 
0.001) and histological type (P = 0.04). 

Prognostic significance of ABCC3 expression. ABCC3 -positive 
expression correlated with shorter OS of NSCLC patients (P log- 
rank < 0.001) (Fig. 4). In univariate analysis, positive pN status, 
advanced TNM stage, and ABCC3 -positive expression correlated 
with poor prognosis (P < 0.001). In multivariate regression model, 
ABCC3 expression remained a significant prognostic factor inde- 
pendent of pN status, pTNM stage and histological type (P < 
0.001) (Table 2). 

Expression of ABCC3 and in vitro drug sensitivity. Among the 199 
samples, 69 (34.6%) were sensitive, 44 (22.1%) were resistant, and 86 
(44.2%) were partially sensitive to the 5 anticancer drugs. 15 (21.7%), 




sensitive group 



resistant group 



Figure 2 | ABCC3 expression by real time-PCR in NSCLC. ABCC3 
mRNA expression was higher in the resistant group compared to the 
sensitive group (P < 0.001). 
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Figure 3 | ABCC3 protein expression in NSCLC. Positive (panel A-B) and negative (panel C-D) expressions of ABCC3 in patients were observed. (A, C, 
X100, SP; B, D, X200, SP). 



20 (23.2%) and 26 (59.1%) samples of the sensitive, partial sensitive 
and resistant group, respectively, were ABCC3 positive. A significant 
association was observed between ABCC3 -positive expression and 
reduced sensitivity to the 5 anticancer drugs using the MTT assay (P 
< 0.001), suggesting that ABCC3 -positive expression is predictive of 
chemoresistance to the 5 anticancer drugs (Supplemental Table 4). 

Discussion 

In this study, we performed in vitro chemosensitivity tests using 
"new generation" anticancer drugs for 199 NSCLC patients. We 



observed 212 genes were significantly altered in the resistant group 
compared to the sensitive group. These genes are involved in various 
biological processes such as cell apoptosis, cell cycle, cell adhesion, 
signal transduction and ABC transporter, which have been shown to 
predict chemosensitivity in NSCLC. The resistant group with 
ABCC3 -positive expression had poor OS. 

Multiple transporters constitute the major mechanisms of 
MDR 8 " 11 . ABC transporters confer resistance to many anticancer 
agents. We identified ABCC3 which has been associated with the 
prognosis of different type of tumors, as one of the most upregulated 



Table 1 Association between ABCC3 expressions an 


d clinicopathologic factors in patients 


with NSCLC 




Variable 


No. of Patients 


(-) 


ABCC3 Expression 
( + ) (++) 


(+ ++) 


Positive Rate (%) 


P(x 2 ) 


All cases 


199 


136 


28 


23 


12 


31.66 




Age (years) 
















<60 


102 


69 


12 


15 


6 


32.36 


0.456 


>60 


97 


67 


16 


8 


6 


30.93 




Gender 
















male 


106 


74 


17 


9 


6 


30.19 


0.459 


female 


93 


62 


1 1 


14 


6 


33.33 




pT status 
















pTl 


78 


57 


1 1 


6 


4 


26.92 


0.146 


P T2 


108 


72 


16 


15 


5 


33.33 




pT3 


13 


7 


1 


2 


3 


46.15 




pN stage 
















pNO 


76 


61 


1 1 


1 


3 


19.74 


<0.001 11 


pNl 


81 


51 


16 


10 


4 


37.04 




pN2 


42 


24 


1 


12 


5 


42.86 




pTNM stage 
















1 


72 


58 


10 


1 


3 


19.45 


<0.001 11 


II 


76 


50 


17 


8 


1 


34.21 




III 


51 


28 


1 


14 


8 


45.10 




Histological classification 
















Squamous cell carcinoma 


91 


68 


5 


12 


6 


25.27 


0.04^ 


Adenocarcinoma 


96 


58 


21 


1 1 


6 


39.59 




others 


12 


10 


2 


0 


0 


16.67 




T P value was from Fisher's Exact test. 
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more malignant phenotype in lung ADC patients. ABCC3 overex- 
pression has been implicated in acquired MDR in cancer cell lines. 
ABCC3 also contributes to resistance to doxorubicin in breast can- 
cer 23 , to 5-fluorouracil in pancreatic cancer cells 24 and resistance to 
adriamycin in hepatocellular carcinoma cells, which however, can be 
effectively reversed by tetramethylpyrazine 25 . Our results are consis- 
tentwith ABCC3 expression being predictive of resistance to pacli- 
taxel, docetaxel, gemcitabine, vinorelbine and cisplatin in NSCLC. 

The limitation in our study includes the fact that the results of 
MTT were not compared with the optimal regimen. The current 
study did not include the functional study of drug- resistant related 
genes that have been screened, thus further studies are needed to 
explore the detailed process of multiple-drug resistance and to pro- 
vide more evidence of how ABCC3 is involved in MDR in vitro. 
Although we collected samples from large amounts of NSCLC 
patients, only 1 1 samples were eventually selected for gene express- 
ion microarray, which also limited the extent of application. Our 
observations need to be validated in a multiple- centered study invol- 
ving patients with a large sample size. However, our observations 
support the notion that chemotherapy combined with targeted ther- 
apy may improve the efficacy in patient with MDR tumors. 

Most previous studies have screened drug- resistant genes from cell 
lines in vitro, while our study validated and predicted drug resistance 
by using primary cells of human cancer tissues. We conducted drug 
resistance experiments, identified drug- resistant genes directly from 
clinical specimens of NSCLC patients, and validated the observations 
in vivo, which accurately reflected the original drug sensitivity vs. 
resistance in patients. We have employed a new method to popular- 
ize the MDR genes screening, which has clinical implications for 
individualized therapy for a variety of tumors. 

In summary, MTT drug response assay in surgically resected 
NSCLC may have clinical utility for predicting tumor response to 
chemotherapy and to prospectively identify drug sensitive and res- 
istant subpopulations. ABCC3 -positive expression correlated with 
enhanced chemoresistance and reduced OS in NSCLC patients. 
For patients with resectable stage IB to IIIA NSCLC, cisplatin-based 
regimen (in combination with Paclitaxel, Docetaxel, Gemcitabine, 
Vinorelbine, etc.) is the majorchoice for postoperative adjuvant ther- 
apy. However, some patients develop multiple-drug resistance, 
which leads to failures in therapy and a poor prognosis. By using 
in vitro MTT assay, we could evaluate an individual patient's sens- 
itivity vs. resistance to conventional anticancer drugs before the 
treatment is employed. Based on the results, an unnecessary cyto- 
toxicity/side effects caused by an insensitive drug could be avoided. 
In our study, we found the correlation of ABCC3 expression with 



Table 2 Prognostic factors 
Variables 


for NSCLC patients 
HR 


Univariate 95% CI 


P 


HR § 


Multivariate 95% CI 


P 


Age (year) 














<60vs. >60 


1.034 


0.584-1.830 


0.909 


_§ 






Gender 














Male vs. female 


1.291 


0.733-2.275 


0.377 








pT status 














pT2-4 vs. pTl 


1.009 


0.621-1.639 


0.971 








pN stage 














pNl-2 vs. pNO 


2.388 


1.640-3.477 


<0.001 








pTNM stage 














ll-lll vs. 1 


2.461 


1.691-3.583 


<0.001 








Histological classification 














Adenocarcinoma vs. Squamous cell 0.653 


0.403-1 .058 


0.083 








cancer 














ABCC3 expression 














ABCC3(+)vs. ABCC3(-) 


2.686 


2.077-3.472 


<0.001 


2.635 


1.992-3.486 


<0.001 


§ HR, 95% CI and Pvalues were from mu 
§ The variable with a P § value of > 0.05 


tivariate Cox proportional ha 
n the multivariate model was 


zard regression model including al 
not shown. 


the demographic and 


clinical factors. 
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Figure 4 | Kaplan-Meier analysis of OS in 199 NSCLC patients according 
to expression of ABCC3. The patients with ABCC3 -positive expression 
had reduced OS. 

genes in the chemoresistant group. In pancreatic cancer ABCC3 
mRNA has been shown to be upregulated and associated with tumor 
grading 17 . ABCC3 overexpression in glioblastoma multiforme 
(GBM) also correlates with a higher risk of death and may be a 
potential therapeutic target for GBM 18 . In breast cancer amplification 
of ABCC3 occurs in primary breast tumors, predominantly in HER2- 
amplified and luminal tumors 19 . In our study, 31.6% of the NSCLC 
samples expressed ABCC3 and its expression correlated with 
advanced pN status and pathological TNM stage and the patients 
with ABCC3 -positive expression had 2.63 times higher risk for death. 

ABCC3 was highly expressed in tumor cells but not in type II 
alveolar pneumocytes 20 . Moreover, the number of ABCC3 -positive 
cells was significantly increased compared to the precursor lesion of 
lung adenocarcinoma (atypical adenomatous hyperplasia, AAH) to 
localized non-mucinous bronchioloalveolar carcinoma (LNMBAC). 
Together, these data suggest that ABCC3 is a candidate novel tumor 
marker for histological diagnosis and functions as a driving factor in 
carcinogenesis. We found ABCC3 expression was higher in ADC 
than SCC. Recent studies have shown that gene expression patterns 
of lung ADC differ substantially from SCC, and each has different 
etiology, carcinogenesis, and clinicopathological behavior 21 ' 22 . These 
observations suggest that high expression of ABCC3 may confer a 



SCIENTIFIC REPORTS | 3 : 3 1 20 | DOI: 1 0.1 038/srep031 20 



4 



multiple-drug resistance and reduced overall survival, confirming 
that in vitro MTT assay is reliable and feasible in NSCLC. 

Therefore, we propose ABCC3 as a novel marker to predict tumor 
response to chemotherapy and postoperative clinical outcome of 
NSCLC patients, which will be helpful to stratify patient subpopula- 
tion for an efficient therapy and realize the concept of personalized 
medicine. Future studies will investigate whether the inhibition of 
ABCC3 increases the sensitivity to chemotherapeutic agents and 
provides new effective targeted therapy of NSCLC. 

Methods 

Patients and tissue specimens. Each patient had signed informed consent for 
medical record review and tissue sample donation before surgery. This study was 
approved by the institutional review board at Harbin Medical University and 
conducted according to all current ethical guidelines. Fresh tumor tissues from 205 
patients with NSCLC were obtained at the Department of Thoracic Surgery of the 
Affiliated Tumor Hospital of Harbin Medical University between May 2008 and 
August 2009. The patients were follow-up until June, 2012. 

Histopathologic subtypes were determined according to the WHO classification. 
The International Staging System was used for pathological staging. All tumors were 
staged according to the pathological tumor/node/metastasis (pTNM) classification of 
the American Joint Committee on Cancer (7th edition). All tumors were histologi- 
cally reviewed by pathologists. Tumor tissue collected at the time of curative surgical 
resection was placed immediately into normal saline. Then a pathologist conducted a 
representative tumor biopsy and stored the tissue in RPMI-1640 for primary culture 
of NSCLC cells and subsequent MTT assay within 30 minutes. Long-term storage of 
the tissues was at — 80°C. At the time of surgery and tissue collection, none of the 
patients had received neoadjuvant chemotherapy or radiation therapy. 

Primary culture of NSCLC cells. Each fresh specimen was homogenized and 
digested in 0.25% trypsin. The dispersed cancer cells were filtered through an 80 um 
nylon mesh, washed twice in 500 uL RNase free PBS, centrifuged at 1000 rpm for 
5 min, and then incubated in RPMI 1640 medium supplemented with 10% fetal 
bovine serum and 1% (v/w) penicillin/streptomycin at 37°C in a humidified 
atmosphere of 5% C0 2 . The samples of 6 patients were contaminated and excluded 
from the consequent experiments. In total, the culture of primary cells was successful 
in 199 samples. 

MTT drug response assay. The NSCLC cells from each patient were seeded at a 
density of 1 X 10 5 per well in 96-wellplates with 190 uL culture medium for 24 h. 
Then the cells were treated with 10 uL gradient concentrations (10X, 5X, 2.5X, 
1.25 X, 0.625 X, 0.3 12 X, 0.156 X ppc) of each of the drugs paclitaxel, docetaxel, 
gemcitabine, vinorelbine and cisplatin for 72 h. "ppc" represents peak plasma 
concentration (the highest level of drug that occurs in the plasma portion of blood 
after administration of the drug). The tumor cells treated with PBS served as controls. 
Then 20 uL tetrazolium salts (MTT, Sigma, 5 mg/mL) were added to each well and 
incubated for 4 h. The salts were cleaved into a colored formazan product by 
metabolically active cells and 100 uL DMSO was added. 

After checking for complete solubilization of the purple formazan crystals, the 
optical density (OD) was measured at 490 nm wavelength. The inhibition rate (IR) 
was calculated as the formula: IR = (OD control group - OD exp erimentai group)/ 
OD control group X 100% 26 " 28 . IR > 50% was defined as highly sensitive; IR 30% ~ 50% 
as moderately sensitive; and IR < 30% as resistant. The data were presented as mean 
± standard deviation (SD) from four replicate wells per microtiter plate. The samples 
with the highest sensitivity or resistance to all the 5 drugs were selected for gene 
expression microarray. 

Gene expression microarray. Extraction and purification of total RNA was 
performed using TRIZOL reagent (Invitrogen, Carlsbad, CA). RNA purity and 
integrity was controlled by Bio Analyser 2100 (Agilent technologies, CA, USA). Gene 
expression profiling was performed using Human Genome U133 Plus2.0 Affymetrix 
(Santa Clara, USA), which included 47,400 transcripts and represented 38,500 human 
genes. Target preparation for Affymetrix Genechips™ was performed according 
tothe manufacturer's instruction. 4 |ig of total RNA was reverse-transcribed and 
amplified by using the RNA amplification kit (Ambion, TX, USA). Hybridization, 
washing, staining and scanning of chips were carried out according to standard 
protocols (Affymetrix, Santa Clara, CA). The mean ± SD of replicated spots was 
calculated for each geneby using the Acuity 4.0 software (Molecular Devices, LLC, 
USA). 

Microarray validation. Real time-PCR. The mRNA level of each gene was measured 
by real-time- PCR using TaqMan method and quantified using an ABI PRISM 7000 
Sequence Detection System (Applied Biosystems, Foster City, USA). GAPDH was 
used as the internal control. All experiments were performed in triplicate. Relative 
mRNA abundance was calculated as 2" Act [Act = Ct (ABCC3) " C « GAPDH >]. 

Immunohistochemistry (IHC). IHC assay was performed using the labeled strepta- 
vidin-biotin (LSAB) method. In brief, tissue sections (4 urn) were cut from the tissue 
blocks. To retrieve antigenicity, tissue slides were immersed in 1 mM EDTA (pH 8.0) 



and bathed in a steamer at 100°C for 15 min. To block the endogenous peroxidase 
activity, the slides were incubated in methanol containing 0.3% H 2 0 2 for 15 min. The 
slides were washed with PBS and incubated with a specific mouse monoclonal 
primary antibody for human ABCC3 ( 1 : 800 dilution, Abeam, Cambridge, UK) at 
4°C overnight. After washing with PBS, the slides were incubated with streptavidin- 
biotinylated anti-mouse IgG (Dako, Carpinteria, CA, USA) as the secondary antibody 
for 10 min at room temperature. The slides were washed with PBS for 3 times and 
stained using diaminoben-zidine tetrahydrochloride (DAB), counterstained using 
hematoxylin, dehydratedand placed on coverslips. The negative controls were per- 
formed by omitting the primary antibody and replaced with PBS. Adenocarcinoma of 
intestine was used as a positive control, which was provided with kit and confirmed to 
express this protein. The positive staining was observed in the membrane of a cell. 

Scoring. The staining results were interpreted independently and blindly by two 
pathologists. For each slide, three to five randomly selected fields were evaluated. For 
each field, the percentage of DAB-positive-tumor cells was calculated as: [(number of 
DAB-positive tumor cells/total number of tumor cells) X 100]. The relative staining 
intensity was defined as negative for <5%; weak ( + ) for 5-25%; moderate ( + + ) for 
25-50%, and strong ( + + + ) for >50% of the tumor cells stained positive for ABCC3. 

Statistical analysis. The association between ABCC3 expression and tumor 
characteristics (tumor size, lymph node, pathological TNM stage, histological 
classification) or sensitivity to chemotherapy was evaluated by % 2 test and Cox 
regression. Kaplan-Meier plot was used to estimate median survival time (MST). 
Overall survival (OS) time was calculated from the date of diagnosis to the date of 
death or last follow-up. The association between ABCC3 expression and survival time 
was analyzed by log- rank test. The effect of ABCC3 expression on OS was estimated 
using hazard ratio (HR) and 95% CI was estimated using Cox proportional hazard 
regression model, with adjustment for other clinical factors. A P-value <0.05 was 
considered statistically significant. All analyses were performed using SPSS 18.0 
(SPSS Inc., Chicago, IL). 
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